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Abstract 

Previous works have been devoted to the study of two-dimensional noninvert- 
ible maps, obtained using a couphng between one-dimensional logistic maps. This 
paper is devoted to the study of a specific one, in order to complete previous re- 
sults [Zj, regarding the evolution of basins and attr actors, when considering the 
tool of critical manifolds. 

1 Introduction 

Many biological and physical phenomena are due to a complex connectivity among 
complex elements; feedback mechanism seems to be a suitable way to understand such 
phenomena (see for instance The introduction of feedback through a control 

parameter leads to models involving global multiplicative coupling among maps. One 
of the most simple models can be given by using logistic maps. Within this context, 
two-dimensional noninvertible maps with a global multiplicative coupling between one- 
dimensional logistic maps are of great interest. The understanding of their quite complex 
behaviour can be a first step in the understanding of dynamics in higher dimensional 
systems. Previous works (see P[S1|E][ZI) have already been devoted to such kind of 
maps. In this paper, our aim is to complete the understanding of dynamics of model (c) 
given in 

The considered model is the noninvertible two-dimensional map T defined by: 

r Xk+i = \{Syk + l)xk{l - Xk) 
\ Vk+i = A(3xfc+i + l)l/fc(l -Vk) 

where A is a real control parameter, x and y are real state variables. The map (0) 
includes a time asymmetric feedback. Numerical studies show that the system ((H) is 
stable when A G [-1.5, 1.19]. 

Section 2 is devoted to recalls about fixed points of the map (^J and their bifur- 
cations. Section 3 concerns critical manifolds. Section 4 analyzes evolution of basins 
by considering crossing of basin boundary through critical curves. In section 5 we ex- 
plain evolution of invariant closed curve issued from a fixed point and giving rise to 
weakly chaotic ring (WCR), then to chaotic attractor. 
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2 Fixed points 

As explained in [B], T possesses at most five different fixed points : 

Po : = 0, 2/0 = (2) 

Pi-xi = " ^2 : a;2 = 0, y2 = (3) 

l-^/a , 4A-3 , 3 

P3-X3 = — - — , 1/3 = — - — , where a = — - — , A > - (4) 

P^:X4 = 2/4 = A > - (5) 

For — 1 < A < 1, Pq is an attractive node and Pi and P2 are saddle points on the 
negative side of the axes, x-axis being Pi unstable manifold and y-axis being P2 unstable 
one. When A = 1, Pq = Pi = P2. For A > 1, Pq is a repulsive node and Pi and P2 
are saddle points, but located on the positive side of the axes, which are now Pi and 
P2 stable manifolds (Figures EHZI)- When A = —1, Pq undergoes a fiip bifurcation and a 
period doubling cascade occurs when A decreases. P3 and P4 exist when A > |. When 
A = |, P3 = P4; for I < A < 1, P3 is a saddle point and P4 an attractive node, the 
whole diagonal segment between P3 and P4 is locus of points belonging to heteroclinic 
trajectories connecting the two fixed points. When A = 1, P4 undergoes a Nei'mark-Hopf 
bifurcation and becomes repulsive by giving rise to an attractive invariant closed curve 
(ICC) (C). Section 5 is devoted to the study of (C) evolution. 



3 Critical lines 

An important tool used to study non-invertible maps is that of critical manifold, which 
has been introduced by Mira in 1964 (see [111101 for more details). A non-invertible 
map is characterized by the fact that a point in the state space can possess different 
number of rank-1 preimages, depending where it is located in the state space. A critical 
curve LC, in the two-dimensional case, is the geometrical locus of points X in the state 
plane having two coincident preimages, T~^(X), located on a curve LC-i. It is recalled 
that the set of points T~"(X) constitutes the rank-n preimages of a given point X. The 
curve LC verifies r~^(LC)^LC_i or T(LC_i)=LC and separates the phase plane in two 
regions where each point has a different number of preimages (k in one region and k+2 
in the other one, k G N). For the map (Q), LC^i is given by the cancellation of the 
Jacobian of the map. 

For the map (^, the critical manifold LC_i is given by four arcs (Figure Q): 

Lc_ia : 2; = ^, Lc_u ■y = ^, Lc_id ■y = -\ (6) 
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^^-^^ '-^ = -1- 9Ax(l-x) 
and the critical manifold LC by (Figure |2)): 

i^a:2/ = y(3x+l)(^-l)(l-^) (8) 
A 

LC,:y=-{?>x+l),x<— (9) 

i^Ce= |(a: = -^,y = 0)| cLCnLCb (10) 

LCd = ((x = 0,2/ = C LCa (11) 



Let us call Qc={x = = ji'^ — h 1)), Qc is a cusp point of LCb, indeed LCb is 

constituted of two merged critical half straight lines {x G ] — oo, -y]), which explains 
the fact that, by crossing LCb, one goes from area to Zq one, or from Z2 to Z2 (cf. 
Figure 121). Let us remark that LCa fl LCb = {Qc, = — |, y = 0)}. 

The LC curves separate the state space in different areas Zj, i = 0,2,4 (Figure 
121), where each point possesses i distinct rank-1 preimages. Figure El shows the state 
space as a foliated one, each sheet corresponding to the existence of a rank-1 preimage. 
This foliation looks like those given for parameter plane, for instance in |H1. Previous 
representations of such foliation related to the existence of LC curves can be found in 

4 Evolution of basins 

The basin of an attractor B is defined as the set of all initial conditions that converge 
towards an attractor at finite distance when the number of iterations of T tends towards 
infinity. The set B verifies : T~^{B) = B and T{B) C B. We call immediate basin the 
largest connected part of a basin containing the attracting set. In this section, we study 
the evolution of basin of attractors when parameter A is modified. The situations are 
shown in Figures El to fTTl Using the terminology and the results of lljjinj; our aim is to 
explain how basins become fractal, nonconnected or multiply connected in the case of 

©• 

• Figure El shows the first modification of the basin with the creation of a bay H, 
due to the crossing of LCa through the basin boundary. The bay H is a rank-1 
preimage of the area which is located inside Z2, when A increases from the 
value 1. 
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• The second modification consists in a bifurcation connected basin ^ nonconnected 
basin wlien A decreases from the value 1. Islands are created (cf. Figures EHZI) after 
the bifurcation concerning the fixed points Pq, -Pi, P2- Indeed, when A decreases 
from the value 1, Pi and P2, located on positive part of axes, go to negative part of 
axes; the stable manifolds of these two saddle points is a part of the basin bound- 
ary; when the two points are on the negative part of axes, a tongue is created in 
the area (x < 0, y < 0) of the state space and crosses through LCa and LCh curves 
to penetrate inside Z2 area, which is located above LC^; it gives rise to a sequence 
of infinitely many preimages of any rank of this part of the tongue: these are the 
islands, the basin becomes non connected. 

• The third modification of the basin is due to a sequence of aggregations of islands 
(Figures IS HTTj) . The first aggregation occurs when LCa and LC^ curves become 
tangent in points Qa and Qb at the basin boundary (Figure ITn|l . rank-1 preimages 
of Qa and Qh-, Qa^ and Q^^, correspond to the aggregation of two islands with the 
immediate basin. The other preimages of Qa and Qb give rise to other aggrega- 
tions. This first aggregation is followed by infinitely many new aggregations due 
to tangencies between critical curves and the immediate basin boundary when A 
decreases; such phenomena are explained more precisely in The important 
point to note is that a crossing between critical curves and basin boundary makes 
some domains change Zi area, so it makes appear or disappear preimages of the 
considered domains. The sequence of aggregations is nearly finished in Figure 
im The basin size increases when the parameter A decreases and tongues corre- 
sponding to island ends disappear when crossing critical curves of different rank 
(Figure IT^ . In Figure IT!^ the basin boundary is smooth. When A decreases, being 
negative, an inverse process gives rise to new appearance of tongues (Figure E}- 

• Another modification corresponds to a bifurcation connected basin ^ multiply — 
connected basin, that means basin with holes inside. The appearance of these holes 
is also related to crossing of basin boundary through critical curves; for instance, 
the hole in the middle of the basin in Figure ^1 is due to the crossing of LCa 
near its left extremum; it makes appear an area in Z2 region and the preimages 
of any rank of this area gives rise to infinitely many holes inside the basin; these 
holes accumulate along the basin boundary. When A decreases, the size of holes 
increases (Figure ITH|). Some holes can open and be changed in bays (cf. [^pil]). 
this is the case in Figure El the basin boundary becomes fractal. 



5 Evolution of attractors 

This section is devoted to the study of evolution of attractors when A increases. This 
evolution can involve relations with the attractor basin or homoclinic and heteroclinic 
bifurcations. 
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• When A= 1, the attractive focus undergoes a Neimark-Hopf bifurcation and be- 
comes repulsive when giving rise to an attractive ICC (C) fFigurelTHj). Oscillations 
appear on (C) when it is tangent to the cusp point Qc on LCb (Figure IT9|l. From 
A = 1.11, A increasing, there is an alternation between existence of ICC with 
oscillations, frequency lockings, cyclic chaotic behaviors, contact bifurcations with 
basin boundaries and single chaotic attractor. These results are consistent with the 
Liapunov exponents values obtained in 0. We are going to explain what occurs 
regarding contact with basins of coexisting period-3 attractors. 

• In the same time, two attractive period-3 orbits appear simultaneously by a fold 
bifurcation, coexisting with two period-3 orbits of saddle points (Figure UHl). So, 
three different attractors coexist; Figure QUI shows the basin of the ICC and the 
basin of the two period-3 attractors. When A increases, the two period-3 attractive 
cycles undergo a Neimark-Hopf bifurcation and become unstable by giving rise to 
two period-3 attractive ICC around them. Figure |^ shows the three different 
basins of the three coexisting attractors : in this case, two period-3 ICC and a 
period-29 attractive cycle, which corresponds to a frequency locking for the ICC. 

• When A increases, the ICC appears again; curies arise, due to crossings of (C) 
through critical curves, and involve auto-intersections of (C), the ICC is then 
called a weakly chaotic ring (WCR) (TU] (Figure |211). The existence of curies 
changes the shape of the ICC and generates contact with its basin boundary; this 
implies the explosion of (C) and makes it disappear. Then, the basin of the former 
ICC becomes part of the basin of the other attractors (Figures 1221 fc . 

• Figure 1221 shows the separate basins of each part of the two period-3 attractive 
ICC, that means basins for are plotted. 

• Another contact bifurcation occurs, which involves the manifolds of the period-3 
saddle cycles, born by fold bifurcation. This contact creates a new attractor, by 
mixing the two former period-3 ICC (Figure EH)). This new attractor evolves, 
undergoes frequency lockings, which give rise to the appearance of cyclic chaotic 
attractor (Figure I77jl : curies appear on the attractor and it becomes a WCR 
(Figures |2HH^); then it becomes a chaotic attractor; chaos becomes stronger and 
stronger (Figures E2HSS1)- 

• The chaotic attractor disappears when a contact bifurcation occurs between the 
attractor and its basin boundary (Figure I^J. 

6 Conclusion 

Complex behaviour of a two-dimensional noninvertible map with a multiplicative global 
coupling between one- dimensional logistic maps, including a time asymmetric feedback. 
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has been analyzed. The different kinds of basin bifurcations which occur when the 
parameter is modified and the evolution of attractors issued from Neimark-Hopf bi- 
furcation have been explained. The fundamental point, as previously obtained in two- 
dimensional noninvertible maps studies, is the position of basins and attractors, regard- 
ing critical curves. Islands aggregations, connected basin nonconnected basin, 
connected basin multiply — connected basin, change of Invariant Closed Curve into 
Weakly Chaotic Ring or Strong Chaotic Attr actor are bifurcations in relation with 
critical curves. 
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lambda = 1 



Figure 1: Preimages of critical lines: 
LC_i arcs. 



Figure 2: Critical lines LC and Zi areas in 
the state plane. 




Figure 3: Qualitative representation of the foliation of the state space, the open U- 
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lambda = 1 
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Figure 4: A=l, the basin is a square. 



Figure 5: Creation of a bay H when LCa 
crosses through basin boundary; A>1 




Figure 6: Creation of islands. 
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Figure 7: Before the bifurcation 
connected basin ^ nonconnected basin. 



lambda = 0.72 
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Figure 8: Before islands aggregation. 




Figure 9: Islands aggregation. 



Figure 10: Enlargment of Figure IHl 
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lambda = 0.2 




Figure 11: A sequence of islands 
aggregations occurs when A decreases. 



Figure 12: Tongues disappear by crossing of 
basin boundary through critical curves. 
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i7- 1 o I ■ I J • j-u Figure 14: Tongues appear again when A<0. 

ligure 16: Ihe basm boundary is smooth. * o i-^- o 
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lambda = -1.1 




Figure 16: When A decreases, the size of 
Figure 15: Basin becomes multiply holes increases. Pq has undergone a flip bi- 

connected. furcation: an attractive period-2 cycle ex- 

ists. 




Figure 17: Basin boundary is fractal. An attractive period-2 cycle exists. 
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Figure 18: An attractive invariant closed 
curve (C) exists around the fixed point P4. 
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Figure 19: Two stable 3-periodic orbits ap- 
pear by fold bifurcation; oscillations appear 
on (C) when it is tangent to the cusp point 
Qc on LCfo. 
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Figure 20: Basins of the ICC and of the 
two period-3 attractive cycles. 



Figure 21: Frequency locking on a period-29 
cycle for the ICC and Neimark-Hopf bifur- 
cation for the two period-3 attractive cycles. 



Basin bifurcations in a 2-D map 



15 



lsmbds= 1.1172 



lambda = 1.12959 




Figure 22: Transient after that the ICC 

(C) has disappeared, the initial condition Figure 23: Basins for the map T^. It per- 

(.5, .51) gives rise to a trajectory, which mits to see the basin of each curve, part of 

follows the shape of the former ICC (C), period-3 ICC. 
before reaching one of the period-3 ICC. 



l2ir*da = 1.117 IsiT^rlz = 1.1 17 




Fiffure 24: The curve (C) becomes a ine.akhi r.hantir. rinn when curies arise on it. 
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Figure 25: Two 3-periodic stable ICC Figure 26: After contact bifurcation, a sin- 
and their basin. gle attractor exists. 
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Figure 27: Cyclic chaotic attractor. 



Figure 28: Oscillations increase. 
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Figure 29: Curies issued from heteroclinic 
tangencies. 
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Figure 30: Enlargment of Figure 
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Figure 31: Oscillations increase on the 
WCR. 



Figure 32: Chaotic attractor. 
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lambda = 1.19 




Figure 34: Contact bifurcation of the 
Figure 33: Chaotic attractor. chaotic attractor with its basin boundary. 



